Astronomy <fe Astrophysics manuscript no. H4220 


February 2, 2008 


(DOI: will be inserted by hand later) 





Amplification of MHD waves in swirling astrophysical flows 

Andria D. Rogava 1 ' 2 *, Gianluigi Bodo 3 , Silvano Massaglia 1 , and Zaza Osmanov 4 



O 
O 



(N 



> 

o 

m 
o 

o 

Or 

6 ■ 

a: 

<ya . 



1 Dipartimento di Fisica Generale, Universita degli Studi di Forino, Via Pietro Giuria 1, Forino F10125, Italy 

2 Abdus Salam International Centre for Fheoretical Physics, Frieste 1-34014, Italy 

3 Osservatorio Astronomico di Forino, Strada dell'Osservatorio 20, 1-10025, Pino Forinese, Italy 

4 Centre for Plasma Astrophysics, Abastumani Astrophysical Observatory, Fbilisi-380060, Georgia 



Received / Accepted 



Abstract. Recently it was found that helical magnetized flows efficiently amplify Alfven waves ( |Rogava et al. 2003 1 . 
Fhis robust and manifold nonmodal effect was found to involve regimes of transient algebraic growth (for purely 
ejectional flows), and exponential instabilities of both usual and parametric nature. However the study was made 
in the incompressible limit and an important question remained open - whether this amplification is inherent to 
swirling MHD flows per se and what is the degree of its dependence on the incompressibility condition. In this 
paper, in order to clear up this important question, we consider full compressible spectrum of MHD modes: Alfven 
waves (AW), slow magnetosonic waves (SMW) and fast magnetosonic waves (FMW). We find that helical flows 
inseparably blend these waves with each other and make them unstable, creating the efficient energy transfer 
from the mean flow to the waves. Fhe possible role of these instabilities for the onset of the MHD turbulence, 
self-heating of the flow and the overall dynamics of astrophysical flows are discussed. 



1. Introduction 

Lately it became quite commonly presumed and hoped 
that swirling three-dimensional motion occurs in different 
kinds of space plasma flows. Fhe presumption is grounded 
on the mounting evidence from different branches of ob- 
servational astronomy, while the hope is related to the 
long-standing aspiration of theoreticians to find in astro- 
physical flows natural laboratories for testing of their ex- 
| pectations and ideas. 

One of the most remarkable observational indica- 
tions came from Solar and Heliospheric Observatory — 
Coronal Diagnostic Spectrometer (SOHO-CDS) data 
that led to the identification IjPike fc Mason 1998|) of 
macrospicules having both rotational and jet-like fea- 
tures. Pike and Mason presented evidence for the ex- 
istence of blue- and red-shifted emission on either side 
of an axis of a macrospicule stretching above the limb 
from a footpoint region on the disk. They interpreted 
these observations as indicating the presence of a ro- 
tation within these tall squalls of hot and magnetized, 
swirling plasma flows — christened as solar tornados 
- observed both on the limb and the disk of the 
Sun. This discovery allows to argue (for earlier argu- 
ments in the same vein see IPneuman fc Orrall 19861 
Shibata & Uchid a~P986|l that the wide range of dynamic 
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events in the solar atmosphere, including micro flares, jets, 
plumes, surges, spicules and macrospicules, may exhibit 
complicated, kinematically nontrivial, three-dimensional 
sheared plasma motions. Generally speaking it looks quite 
credible that rotation plays a major role in the dynamics 
of chromospheric and transition region features. It seems 
plausible to admit that the presence of these complicated 
patterns of plasma motion deeply influence the dynamics 
of MHD waves in the solar atmosphere, contribute to the 
coronal heating, and to the acceleration of the solar wind. 

Another message of observational evidence for swirling 
flows recently came from the totally different class of 
astronomical objects: Herbig-Haro (HH) jets. Two dif- 
ferent observations of young stellar objects in HH 212 
IIDavis et al. 2000f) and DG Tau HBacciotti et al. 2002jl . 
contain serious indications for the presence of rotation 
in these jets. These results are in accordance with pre- 
dictions of the popular magnetocentrifugal jet launching 
model ( Blandford fc Payne 19821 ). The possible impact of 
these complicated, helical motions on physical processes 
occurring within these flows is yet to be understood. 

The growing evidence in favour of astrophysical helical 
flows poses a twofold challenge. At the one hand it nat- 
urally solicits for the further theoretical study of general 
aspects — generation, equilibrium, stability and internal 
dynamics — of helical flows in the framework of plasma 
astrophysics. At the other hand, after acquiring better un- 
derstanding of basic physics, it suggests to build concrete 
prototype models, closely adjusted to specific examples of 
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observed swirling flows, which might appear to be use- 
ful for shedding some light on the puzzling observational 
appearance of related astronomical objects. 

The basic theory of shear flows tells us that dynam- 
ics of waves and vortices, sustained by these flows, are 
substantially affected by the differential character of mo- 
tion. On the basic mathematical level these phenom- 
ena, somewhat misleadingly christened as "nonmodal pro- 
cesses", are related to the non-self-adjointness of linear 
dynamics of perturbations in SF ( Tre fethen et al. 19 93). 
The variety of these processes, that could be labeled as 
"shear- induced nonmodal processes" (SINP), is quite well- 
understood for simple, plane parallel flows and thoroughly 
described in the recent literature (for the latest review see, 
e.g.. lBodo et al. 2001l and references therein). It is known 
that SINP lead to the generation of new modes of plasma 
collective behaviour, to new forms of flow- wave, wave- wave 
and vortex-wave interactions provoked and fueled by 'par- 
ent' shear flows. These phenomena, originally disclosed 
in hydrodynamics, take place in various kinds of plasmas 
and might have a number of astrophysical applications, in- 
cluding pulsar magnetospheric plasmas, solar atmospheric 
phenomena and galactic gaseous disk dynamics. The in- 
terest towards possible applications was strengthened by 
the recent numerical evidence that real-space appearance 
of SINP is easily recognizable and robust even in the pres- 
ence of tangible dissipation l|Bodo et al. 2 001 1. 

Real astrophysical shear flows are almost always in- 
volved in motions with complex kinematics/geometry 
and kinematic complexity is known to bring an addi- 
tional variety to SINP. When geometry and kinemat- 
ics of flows are complex a whole bundle of new ef- 
fects arises, like "echoing" (repetitive) transient pulsations 
and different kinds of shear-induced (including paramet- 
ric) instabilities jMahajan fe Rogava 1 999). The study 
of kinematically complex SF was initiated in hydro- 
dynamics ||Lagnado et al. 1984[ ICraik fc Criminale 1 986 
ICriminaTF'fc^D'razm^M)}) and this approach is still 
largely unknown for the plasma and astrophysics commu- 
nity. Therefore, it is an intriguing challenge and a task 
of a big practical importance to study these processes in 
MHD flows of nontrivial geometry and kinematics. 

The first step in this direction was recently made 
< |Rogava et al. 20 03, hereafter referred as Paper I). In 
this study the simplifying assumption of incompress- 
ibility was adopted, which cuts off modes of acoustic 
origin (slow magnetosonic waves, SMW, and fast mag- 
netosonic waves FMW) and allows to concentrate on 
the investigation of the dynamics of Alfven waves 
(AW). The subject of the interaction between flows and 
AW is interesting in a number of astrophysical appli- 
cations jBalbus fc Hawley 1991[ |Tagger et al. 1992 
Tagger fc Pellat 1999| |Ryutova et al. 2001 

^anuerT^Tagger^002l ) . In Paper I it was found 
that helical shear flows are efficient amplifiers of AW. In 
purely ejectional flows (i.e., when no rotation is present) 
AW are amplified transiently via algebraic, shear-induced 
instability. In a swirling flow AW are exponentially 



unstable: depending on the mode of differential motion 
both usual and parametric instabilities appear. 

These results were discussed in the context of their pos- 
sible (observable) manifestations. It was argued that they 
might account for the generation of the large-amplitude 
Alfven waves e.g., within 'tornado- like' patterns existing 
in the solar atmosphere. It was suggested that they could 
lead to the efficient self-heating of flows: the kinetic energy 
of the flow, being extracted by amplified Alfven waves, re- 
turns back to the flow in the form of the thermal energy 
because AW are eventually damped via magnetic diffu- 
sion. Finally, it was argued that these instabilities might 
serve as an initial (linear and nonmodal) phase in the ul- 
timate subcritical transition to MHD Alfvenic turbulence 
in various kinds of astrophysical shear flows. 

It is known that the incompressibility condition, used 
in the Paper I, is quite restrictive. Often the usage of this 
condition introduces its own imprints on the dynamics of 
perturbations and it is not trivial to distinguish the gen- 
uine SINP from the phenomenological effect imposed by 
the incompressibility approximation. Therefore it is quite 
important to consider MHD waves in helical flows without 
the usage of the incompressibility condition and to study 
the SINP for the full spectra of MHD waves containing 
together with the AW also the SMW and the FMW. 

This task is undertaken in this paper. We find that the 
range of processes, sustained by helical flows is extremely 
rich. It encompasses different kinds of wave transforma- 
tions, wave beatings. Most importantly, usual and para- 
metric shear instabilities, which has been found in the 
Paper I, are found to appear for the full spectra of MHD 
waves containing AW, SMW and FMW. We see that flows 
of such a high degree of complexity efficiently intertwine 
all three MHD wave modes and efficiently exchange energy 
with them. The relevance of these results to the physics 
of helical flows and to the understanding of possible ob- 
servational appearance of related astronomical objects are 
pointed out and critically discussed. 

2. Theory 

Our aim is to study linear collective MHD modes in helical 
flows. For this purpose we need to write equations of the 
ideal linearized MHD for the evolution of perturbations 
within the flow. In Paper I we studied only incompressible 
perturbations, while now we consider fully compressible 
case, so our starting equations are [D t =d t + (V-V)]: 

% + ^(V-V)=0, (1) 

D 4 V = --VP-JLx(VxB), ( 2 ) 
p Airp 

D t B = (B-V)V - B(V-V), (3) 

V B = 0. (4) 

Our equilibrium model, used in Paper I, assumes a 
homogeneous MHD plasma (po = const), embedded in 
a homogeneous, vertical magnetic field (B o =[0, 0, B = 
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const]). We consider instantaneous values of all physical 
variables as sums of their mean (equilibrium) and pertur- 
bational components: B=B +B', p=p +p', etc. Applying 
this decomposition we convert (1-4) into the following set 



for perturbation variables [D t =d t + (U -V)]: 

£> t d + V-u = 0, (5) 

V t u + (u-V)Uo = -C 2 s Vd + C\[d z h - Vb z ], (6) 

D t b=(b.V)U +a,u + e,(V-u), (7) 

Vb = 0, (8) 



with d=p' I po and b=B'/_Bo- Note also that for compress- 
ible perturbations p' — C 2 p' with C s the homogeneous 
speed of sound. 

If we introduce a new vector quantity h=b — e z d then 
the above system can effectively be reduced to the follow- 
ing set of second-order equations: 

V 2 h - [(h-V)U -V]U - {C 2 + Ci)Ah = 

C 2 A [d 2 z h - V(d z h z ) - e z d z (S7-h)}, (9) 

which describes evolution of all three MHD modes - 
SMW, AW and FMW — influenced by the presence of 
the equilibrium flow with an arbitrary Uo- 

In paper I we worked with the velocity field specified 

by 

U(r)=[0, rO(r), U(r)], (10) 

with f2(r) = A/r n , where r = (x 2 + y 2 ) 1 ^ 2 is a distance 
from the rotation axis, while A and n are some constants. 
In particular, n = and A — f^o for rigidly rotating plas- 
mas, while n — 3/2 and A = (GM) 1 / 2 , for the Keplerian 
rotation. 

The linear shear matrix for this kind of U(r) has the 
following form ( |Rogava et al. 2003| ): 

(a Ai 0\ 
S=\ A 2 -a . (11) 
\Ci C 2 Oj 

The key equation of the "nonmodal" method 
IjCraik fc Criminate 19861 |Mahaj an fc Rogava 1999Q , 

d t k + S T k = 0, (12a) 

giving a full evolutionary picture of k's, transcribes in this 
case to the following set of equations 1 : 

4 X) +°k x + A 2 k y + dk z = 0, (126) 

+ A x k x -ak y + C 2 k z = 0, (12c) 

while k z = const. 

These equations imply that , while k x (t) and k y (t) may 
have algebraic, exponential or periodic time dependence. 
For swirling flows the differential rotation parameter n 
plays decisive role in determining the evolution scenario 
for the wave number vector |k(t)|: when n < 1 (including 

1 Hereafter F^ n ' will denote n-th order time derivative of a 
function F. 
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the rigid rotation case) the time evolution of the |k(f)| 
is periodic, while when n > 1 (including the Keplerian 
rotation regime) |k(t)| evolves exponentially. 

Note that the following nonlinear combination of k x (t) 
and k y (t): 

A = k x k y ^ - kykW + k z (dk y - C 2 k x ) = 

A 2 K 2 -A 1 kl + 2ak x k y + 2k z (C 1 ky-C 2 k x ) = const. (13) 

is a conserved quantity. In the case of pure, two- 
dimensional rotation (k z = C\ = C 2 = 0) it reduces to the 
constant A r =fc. c /4 1 ' ) — kyk^ ( |Mahajan fc Rogava 1999| ) , 
while for the case of a pure outflow (A\ = A 2 = a = 0) it 
reduces to the conservation of the quantity A e =Cik y — 
C 2 k x ( Rogava et al. 2 000). Therefore, we see that the 
wave-number invariant A in the case of the helical flow 
is the sum of the A r and the A e functions. 

Using the same method as in Paper I (see for details 
also |M aha j an fc Rogava 1999[ ) we can effectively convert 
the system (5-8) to the set of first oder ordinary differen- 
tial equations [H=ih, D=id\. 

D {1) = k-u, (14) 
U W + Su = -(C 2 + C^)kD+ 

C 2 A [k z tt-kH z +e z k z D], (15) 

H (1) =S-H-/^u, (16) 

k-H = -k z D, (17) 

Note that (14-17) contain two first order ordinary 
differential equations with time-dependent coefficients. 
Time-dependence of these coefficients is completely de- 
termined by the temporal evolution of the wave number 
vector k(i) and, therefore, governed by Eq. (12). We can 
reduce this set to the set of second order equations for the 
components of the vector H, which is the nonmodal form 
of the Eq. (9). It can be written in the following (vector) 
form: 

H (2) + C 2 k(k-H) - S 2 U + C 2 A [(k - k z e z )(k-H) + 

(k 2 z H-kk z H z )} = 0. (18) 

Note that in the helical flow, specified by the shear 
matrix (11), the square of the matrix is equal to 
dRogava et al. 2003| ): 

(T 2 0\ 

||5 2 || = o -r 2 , (19) 
V e x e 2 / 

where we use notation: T=(a 2 +AiA 2 ) x ^ 2 , e\=A 2 C 2 + aC\ 
and e 2 =Aid - aC 2 . Note that ||5 3 || = r 2 ||5||. 

Splitting both H and k vectors into their longitudinal 
and transverse components — H=(H^, H z ),k=(k±, k z ) 
— we can write more explicit form of (18) revealing the 
nature of coupling between the MHD modes, imposed by 
the presence of the velocity shear: 

HP+C^k 2 z H z +C%(k x -U ± )-(e 1 H x +e 2 H y ) = 0, (20) 



4 



Rogava et al.: Amplification of MHD waves in swirling astrophysical flows 



Hi 2) + \C\k\ - T 2 ]H ± + (C 2 S + C 2 A )k ± (k ± -H ± ) + 

C 2 k ± k z H z = 0. (21) 

It is easy to see that in the absence of the shear 
flow Eqs. (20-21) give standard expressions for the disper- 
sion properties of all these modes. The much simpler ver- 
sion of this system was analyzed before for plane-parallel 
SF of standard MHD plasmas ( |Chagelishvili et al. 1996| ). 
In the astrophysical context the similar kind of coupled 
ODE's were studied for MHD waves in the solar wind 
plasmas l|Poedts et al. 1998jl . in galactic gaseous discs 
( |Rogav a et al. 199 9[) and in cylindrical rotationless flux 
tubes jftogava^t a" 2000] ). The factual 'route' of the wave 
number vector k(t) temporal "drift" (caused, in its turn, 
by the differential character of the plasma motion) plays 
the crucial role in the time evolution of physical pertur- 
bations. In parallel flows kj_(t) exhibits linear time de- 
pendence. However, even in this relatively simple case we 
have the whole set of SINP < |Rogava et al. 2000| ): waves ex- 
change energy with the flow and the velocity shear makes 
waves coupled with one another, making possible their re- 
ciprocal transformations. 

Helical flows are expected to exhibit evolutionary 
regimes characteristic to both purely ejectional and purely 
rotational flow patterns. However, we expect that shear- 
induced wave transformations (SIT) are more distinctly 
exhibited by parallel flows without rotation in the trans- 
verse plane, because these flows are stable and SIT is 
the only major kind of SINP occurring in the flow. 
In helical flows, presumably featuring different sorts of 
shear- induced instabilities, SIT could be less- or not- 
pronounced. It should be noted that the parameters e± and 
s 2 exist only when both rotation and outflow are present 
and they are nonzero only when the forces that determine 
the kinematic portrait of the flow are non-conservative 
IjCraik fc Criminate 1986(1 . Thus their role in the dynam- 
ics of MHD waves sustained by swirling flows can be quite 
significant. 

3. Discussion 

The evolution of linear MHD modes in kinematically com- 
plex helical flows might be quite complicated. Before con- 
sidering any special and/or particular cases we have to re- 
call the following three levels of complexity arising within 
this problem. 

First, the medium itself is complex enough, because it 
sustains three different linear modes of oscillations: SMW, 
AW, and FMW. Although in the absence of shear these 
modes are decoupled, still an arbitrary perturbation ex- 
cited within this medium is normally a superposition of 
these three normal modes l|Sturrock 1994). 

Second, in the presence of a simple, plane-parallel flow 
with a linear velocity profile: (a) one mode - FMW - be- 
comes able to draw energy out of the background flow; (b) 
depending on the plasma-/? different waves become cou- 
pled and are able to transform into each other; (c) the sys- 
tem starts exhibiting beat wave phenomena. The resulting 



picture of the MHD wave dynamics becomes considerably 
complex fRogava et al. 2000| ). 

Third, when incompressibility condition is used and 
the presence of the rotation and stretching of flow lines 
in the transverse cross section of the flow is "allowed" 
one finds that new kinds of shear instabilities emerge 
(Ro gava et al. 2003| ). This is expected to happen also for 
compressible (acoustic) wave modes, because even in the 
simplest two-dimensional, nonmagnetized flow with kine- 
matic complexity (Mahaj an fc Rogava 1999| ) usual and 
parametric instabilities do appear together with asymp- 
totically persistent, "echoing" solutions. 

Therefore it seems reasonable to suppose that helical 
MHD flows, possessing all these levels (or degrees) of com- 
plexity plus the complexity of the specifically helical na- 
ture, related with the existence of the E\ and e 2 "heli- 
cal" parameters, must exhibit highly complicated collec- 
tive processes, dominated by different regimes of shear- 
induced variability and instability. 

For numerical purposes it is more suitable to deal with 
the dimensionless version of the (14-17) set: 





= K X {t)V X + K-y(T)Vy 4 




(22) 




-1- T,v x + aiVy = -JC x (t)£ 2 q + b x - 


JC x (r)b z , 


(23a) 




+ a 2 v x — Ev y = —IC y (T)e 2 g + b y - 


JC y (T)b z , 


(236) 




v { p + R x v x + R 2 v y = -e 2 




(23c) 




b x x) = Y,b x + aiby - v x , 




(24a) 




by = a 2 b x - T,b y - v y , 




(246) 




K-x{r)b x + K y (r)by +b z = 


0, 


(25) 



derived from (14-17) with the usage of dimensionless no- 
tation: T=C A k z t, K x (T)=k x (t)/k Zl K,y{T)=ky{i)/k z , 
Vi=Ui/C A , S = <r/C A k z , ai. 2 =Ai. 2 /C A k z , 
Ri^(C h2 /k z C A ), e=C s /C A . 

It is instructive to calculate how the total energy of 
perturbations being, in this case, the sum of kinetic, com- 
pressional and magnetic energies: 

E tot =Ek in + E c + E m , (26) 

E kia =(v 2 x +vl + v 2 z )/2, (27a) 

E c =e 2 g 2 /2, (27b) 

E m =(b 2 x + b 2 y + b 2 z )/2, (27c) 

changes in time. We can easily see that the evolutionary 
equation for the total energy is: 

= (at + a 2 )(b x b y - v x v y ) + X[(v 2 y - v 2 x ) + (b 2 x - 6 2 y )] + 

Ri(b x b z - v x v z ) + R 2 (b y b z - v y v z ), (28) 

The scope of this paper is not full investigation of all 
possible regimes of evolution. Instead, our purpose is to see 
whether shear instabilities, disclosed in the incompressibil- 
ity limit for AW, appear also in the compressible case, for 
the blend of SMW, AW, and FMW. It is instructive and 
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Fig. 1. The temporal evolution of the density D (r) 
and the normalized total energy of perturbations 
£'tot(''")/£'tot(0), which exhibits quasiperiodic and pulsa- 
tional behavior. The set of parameters is: e = 10, /C x (0) = 
10, Ky(0) = 8, Ri = 0.8, R 2 = 2, £ = 0, a x = -0.1, 
a 2 = 0.1. 

convenient to unfold these phenomena by following the 
scheme used in our previous studies of cylindrical flows 
( |Rogava et al. 200tj| and the Paper I) and to see how two, 
major regimes of k(i)-dynamics (periodic and exponen- 
tial) affect the course of collective phenomena in flows 
with different values of the plasma-/?. 

3.1. High-fi plasmas 

In parallel flows, when e 2 ;§>l, FMW are decoupled from 
the AW and SMW, while the latter two are coupled 
< |Rogava et al. 20 00) and may transform into each other. 
Swirling flows show more complicated behaviour. When 
7 2 =(r ICa^z) 2 < 0, i.e., when temporal evolution of 
wavenumber vectors k(t) is periodic, we encounter with 
two new effects: 

1. The appearance of "echoing" waves, consisting of 
repetitive, modulated bumps of AW and SMW, ex- 
changing energy with the mean flow. Fig.l displays 
an example of such process. From the figure we can 
surmise that mutual AW^SMW transformations still 
happen and the resulting mixture of waves "pulsates" 
taking and giving energy from/to the background he- 
lical flow! This process has well-pronounced quasiperi- 
odic nature. 

2. Since coefficients in (22-25) vary periodically it is 
plausible to expect that for certain cases the sys- 
tem must also exhibit some kind of self-parametric 
instability. The term "self-parametric" seems appro- 
priate, because it is the consequence of the veloc- 
ity shear inherent to the system and forcing on itself 
< |Argentina et al.~ 999). This kind of instability was 
first discovered for plain acoustic waves in 2-D flow 



Fig. 2. The evolution of parametrically unstable blend of 
AW and SMW. The set of parameters is the same as on 
Fig.l, except Ri = 0.4. 



Density 




Total energy 




Fig. 3. Shear instability of low-frequency MHD oscilla- 
tions. The plot features temporal evolution of the same 
kind of initial value problem as on Fig.l. The only dif- 
ference is that The evolution of parametrically unstable 
blend of AW and SMW. The set of parameters is the same 
as on Fig.l, except that a\ = 0.1 now, making the sign of 
7 2 positive. 



patterns of neutral fluids ( |Mahajan fc Rogava 1 999). 
Similar sort of instability was found for Alfven waves 
in the Paper I. Numerical examination of the (22-25) 
allowed us to find similar instabilities for the compress- 
ible case. One example is given on the Fig. 2. The values 
of all parameters (except R± ) here are the same as on 
the Fig.l, but R\ = 0.4. The figure shows that in this 
case hydromagnetic oscillations amplify exponentially, 
extracting energy from the flow. Note that for the ex- 
istence of this kind of parametric instability it is not 
necessary to have any periodicity in the background 
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Fig. 4. The temporal evolution of initially excited AW, 
which became partially transformed into FMW and ex- 
hibits quasiperiodic, pulsational behavior in \ow-(3 plasma 
flow. The set of parameters is: e = 0.1, a\ — —0.1, 
a 2 = 0.1, £ = 0, Ri = 0.01 and R 2 = 0.9. Fig. 4 shows 
time evolution of the velocity perturbation v x and the to- 
tal energy normalized on its initial value E io t(r) / '£'tot(0) . 

flow, but it is essential to have periodic time variation 
of the wave number vector. 

Let us turn our attention, now, to flows with "sharper" 
rate of differential rotation, with n > 1. In this case 
7 2 > and the temporal evolution of k(t)'s becomes 
exponential. In the absence of dissipation these flows 
host rather robust shear instabilities. The similar kind 
of instability was found for acoustic waves in a 2-D flow 
(Mahajan & Rogava 1999) and in the Paper I it was de- 
tected for Alfven waves as well. The example of such in- 
stability for the mixture of AW and SMW is shown on the 
Fig. 3. The parameters are the same as for the Fig.l, only 
ai = 0.1 (it reverses the sign of the A 2 ). We see that in 
this case the mixture of AW and SMW undergoes rather 
strong, exponential enhancement of its amplitude and en- 
ergy 

3.2. Low- [3 plasmas 

In this case hydromagnetic oscillations offer quite a dif- 
ferent picture. Namely, the SMW mode is decoupled from 
the other two MHD modes — its dispersion curve runs well 
below the Alfven dispersion "horizontal" and it is not cou- 
pled with other two (AW and FMW) modes. The FMW, 
on the other hand is coupled with the AW. When rotation 
is absent and there is only parallel outflow the coupling 
ensures linear transformation of the AW into the FMW 
(see as an example Fig.6a-d in Rog ava et al. 200"0| ). 

In the helical flow, getting certain "input" of initial 
AW and/or FMW oscillations, these waves will keep trans- 
forming into each other. Besides the mixture of waves 
might exhibit the same kind of "echoing" and unstable 



Fig. 5. The parametrically unstable solution. The values 
of the parameters are: e = 0.1, a± — —0.1 a 2 — 0.1, S = 0, 
Ri = 0.01 and R 2 = 0.83. 
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Fig. 6. The shear instability in the j 2 > case. The values 
of the parameters are: e = 0.1, a± — 0.1 a 2 — 0.1, S = 0, 
Ri = 1 and R 2 = 1. 

behavior as it was seen in high-/? plasmas. Numerical sim- 
ulations support this expectation. Making values of a\ and 
a 2 nonzero (ax = —0.1 and a 2 = 0.1, so that 7 2 < 0) 
and using following parameters: R\ = 0.01 and R 2 = 0.9 
we see (Fig. 4) the appearance of interesting patterns of 
"echoing" solutions with quasiperiodic variability of the 
perturbation total energy 2 . 

We find that parametric instabilities are also charac- 
teristic to low-/3 plasmas. The example is given on Fig. 5. 
The parametric nature of the instability is apparent from 
the remarkable narrowness of the range of parameter val- 
ues, where the instability is present. Note that the set of 
parameters used for Fig. 5 is the same as for Fig. 4 except 

2 Note that here, as elsewhere on the plots of this paper, the 
total energy is normalized on its initial value. 
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Fig. 7. The temporal evolution of the wave blend that 
exhibits quasiperiodic, pulsating behavior in e — 1 plasma 
flows. The set of parameters is: a\ = —0.1 a 2 = 0.1, S = 0, 
i?i = 0.1 and R 2 = 0.85. The plots show time evolution of 
the magnetic field B x perturbations and the normalized 
energy E tot {r) / E tot {Q) . 
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Fig. 8. The evolution of parametrically unstable blend of 
waves. The set of parameters is: a\ = —0.1 a 2 =0.1, £ = 
0, Ri = 0.1 and R 2 = 0.8005. Note that the range of R 2 
when the evolution of the wave mixture is parametrically 
unstable is very narrow: [0.8001; 0.8011]. 



R-2 = 0.83. If one takes the value of R 2 less by 0.01, then 
the instability disappears and the system again displays 
the "echoing behaviour. 

When the differential rotation is characterized by n > 
1 profile (7 2 > 0) the system shows strong exponential 
shear instability. The corresponding plots are given on 
Fig.6. 
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Fig. 9. The shear instability in the A 2 > case. The set 
of parameters is: ai = 0.1 ai = 0.1, S = 0, i?i = 1 and 
R 2 = 1. 

3.3. The case of (3 ~ 1 

From the studies of parallel flows we know that in terms 
of wave couplings and mutual transformations this is the 
most complex case: all MHD wave modes are coupled and 
may transform into each other. In helical flows the pres- 
ence of wave transformations may be less visible, overshad- 
owed (when 7 2 < 0) by quasiperiodic modulation of waves 
appearing as repetitive, ("echoing") bundles of mixed AW, 
SMW and FMW modes (see Fig. 7); or by self-parametric 
instability (see Fig. 8). While when 7 2 > the waves are 
unstable in the similar way (see Fig. 9) as in previously 
considered high-/? and low-/? cases. 

4. Conclusion 

The goal of this work was to find out whether the exotic 
SINP originally found for 2-D kinematically complex ve- 
locity patterns of neutral fluids ( Mah ajan fc Rogava 1999| 
and for 3-D helical flows of magnetized conducting fluids 
in the incompressible limit, appear also for the full spectra 
of MHD waves sustained by compressible MHD medium. 
We found that this is indeed the case! Therefore we can 
now firmly claim that the range of the SINP, typical for 
the flows of complicated, helical nature, is broad. These 
SINP persist to show up both in the incompressible and 
compressible cases. They are present in flows with arbi- 
trary values of plasma-/?. 

In relatively mildly sheared flows (with n < 1, i.e., in- 
cluding rigidly rotating systems) MHD modes appear to 
be rather stable exhibiting either "echoing" pulsational 
behaviour or relatively long-time-scaled parametric insta- 
bilities. One can expect that in helical flows of this nature, 
especially in well-beamed, or well-collimated flows such as 
jets, shear flow effects are not likely to lead to disruptive 
instabilities. Instead, through quasiperiodic interchange of 
energy with the mean flow, they might tend to exhibit cer- 
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tain modes of quasi stable and quasiperiodic structuring 
both in the space and in time. 

In more strongly sheared systems (n > 1 including 
Keplerian rotation) waves within swirling flows become 
subject to potentially very fast-growing shear instabilities, 
which would either lead to the disruption of "parent" flow 
patterns or to the development of the MHD turbulence 
with subsequent phase transition to turbulent rotational 
flow systems. We have to bear in mind that in this case 
the exponential growth of the |k(t)| inevitably makes the 
spatial length-scales of perturbations smaller and smaller. 
It implies that the effects of the viscous decay and/or mag- 
netic diffusion, neglected while we consider the MHD flow 
as an "ideal" one, must sooner or later become important 
and lead to the dissipation of the energy gained by the 
exponentially increasing waves into the heat. It can be 
argued that in accretion-ejection systems, where the rota- 
tional law seems to be quasikeplerian, these instabilities 
may account for the transition to turbulence in accreted 
plasma flows. Alternatively, this process might lead to ef- 
fective "self- heating" of these flows, when the energy ac- 
quired by waves from the flows through the agency of the 
shear instability would eventually transform into heat via 
diffusion. 

Speaking about swirling astrophysical flows we are 
keen to use the term "cosmic tornado" for this class 
of flows, because they are reminiscent of powerful and 
dangerous tornados in the Earth's atmosphere. Recently 
such structures, called solar tornados, were identified 
in the polar regions of the solar atmosphere (both on 
the limb and the disk) by SOHO-CDS observations 
(Pike & Mason 1998). Another class of cosmic tornados 
are, probably, stellar jets, because recent observa- 
tional results (JBacciotti et al. 20021 IDavis et al. 2000|) 
seem to confirm the predictions of various magneto- 
centrifugal jet acceleration models l|Shu et al. 1 995 

IFerreira & Pellet ie r 19951 ICarnenzinri 19~97l 

|Lery et al. 19991 |Konigl fc Pudritz 2000| ) about the 
presence of a swirling motion within the jets. The third 
class of astrophysical flows with possible presence of 
tornado-like motion are accretion columns - magnet- 
ically channeled shear flows of plasma to a neutron 
star's (or a white dwarf's) magnetic pole. The infalling 
matter is decelerated approaching the star surface, but 
it can also form a shock high above the star's surface 
( |Hujeirat fc Papaloizou 1998| ). Originally the formation 
of rotational accretion columns was considered in the 
astrophysical fluid dynamics context fCass en 1978|) . while 
now accretion columns associated with X-ray pulsars and 
cataclysmic binaries are modeled either as thin slabs or 
tall columns of infalling matter. It seems reasonable to 
surmise that accretion columns comprise 3-D swirling 
plasma flows. 

Bearing in mind these perspectives we should stress 
that the results of this paper are quite general and (delib- 
erately) not adjusted to either of three kinds of swirling 
astrophysical flows. For all these classes of 'cosmic torna- 
dos' the level of our factual knowledge about their basic 



kinematic features is still inadequate for building of any 
credible concrete models. However we hope that these fu- 
ture models, based on the gained data and implying real- 
space sophisticated simulations, will show how the generic 
processes disclosed in this paper and the Paper I might in- 
fluence the overall dynamics of real 'cosmic tornados'. 
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